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Abstract 



Quantum mechanical potentials can be divided according to their solv- 
ability into three categories. The best known category is exactly solv- 
able potentials (e.g. the harmonic oscillator and the hydrogen atom), 
which can be solved analytically. This means that all the energy lev- 
els and all the wave functions can be calculated in a finite number of 
algebraic steps. The second category is non-solvable potentials (e.g. 
the quartic potential), for which it is only possible to approximate the 
solution using various analytical and numerical methods. The third 
and least known category is quasiexactly solvable (QES) potentials, for 
which only part of the energy spectrum and the corresponding wave 
functions can be analytically calculated. 

The feature distinguishing QES potentials from non solvable poten- 
tials is that their Hamiltonian may be represented as a block-diagonal 
matrix with at least one finite block. Such a finite block may always 
be diagonalized, thus the energy levels and wave functions correspond- 
ing to it can always be calculated. Block diagonal Hamiltonians can 
be created using bilinear combinations of Lie group generators. The 
exactly solvable potentials may be treated as a sub-group of the QES 
potentials since it is always possible to diagonalizc their Hamiltonian. 

This work deals with the solution of one dimensional QES poten- 
tials. It also includes one dimensional supersymmetric QES potentials 
which are potentials that can be represented by a 2 x 2 matrix and 
their Hamiltonian operates on a wave function which is a two compo- 
nent spinor of one spatial variable. 

First, expanding the results of existing works, some properties of 
the exact solutions of QES potentials are shown to hold for the general 
case. One interesting property of some QES potentials (which is used 
in the second part of this work) is the simple connection between the 
wave functions and the initial value solution of the Schrodinger equa- 
tion. The initial value solution is the generating function of the secular 
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polynomials, the zeros of which are the quasiexactly solvable part of the 
energy spectrum for the potential. This property is shown to hold for 
a class of QES and exactly solvable potentials and not only for specific 
examples. 

Next, the secular polynomials and the QES part of the energy spec- 
trum are approximated using a new technique which is based on the 
WKB approximation. The approximation is useful because calculat- 
ing the energy levels exactly becomes more and more complicated as 
the size of the finite block in the Hamiltonian grows. The initial value 
solution of the Schrodinger equation is calculated using the WKB ap- 
proximation. Then approximate secular polynomials and energy levels 
are derived from the WKB solutions. In order to simplify this process 
saddle point approximation is used to obtain the secular polynomials 
and the energy levels from the WKB approximation. This calculation 
has to be carried out separately for every potential. The process is 
demonstrated for the harmonic oscillator and the sextic potential and 
the quality of the results is discussed. The combined WKB and saddle 
point approximations give a simple evaluation of the energy levels of 
QES potentials with a generating function that solves the Schrodinger 
equation. 
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List of Symbols 



symbol 


meaning 


a{x) 


imaginary phase and exponential behavior of the wave function 


A{x) 


derivative of a{x) and quasigauge field 


Ck 


wave function coefficient 




coefficients of the generators in Hq 


Dk 


even wave function coefficient 


dk 


odd wave function coefficient 


dk 


lower wave function coefficient 


E 


energy 


En 


energy eigenvalue 




highest energy of the QES sector 


H 


Hamiltonian 


Hg 


quasigauge Hamiltonian 


3 


the group parameter 


J 


u{l) group generator 


M2,+l 


coefficient matrix 




minor of the term in the n + 1 row and the k + 1 column of M 




minor of coefficient matrix 




minor of coefficient matrix 


Pk{E) 


secular polynomials 


PiiO 


the coefficients of the ^ derivatives in Hq 


Qk{E) 


even secular polynomials 


qk{E)_ 


odd secular polynomials 


Qa-i Qa 


group generators 


Rj 


the jth representation of a group 


Sa 


phase of Cauchy integral 


S 


scaled x 
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symbol meaning 



S/3 saddle point of 
group generator 

t scaled y 

Uk lower wave function coefficient 

V{x) potential 

v{s) scaled V{x) 

ctfe a coefficient in the recursion of the wave function coefficients 

f3k a coefficient in the recursion of the wave function coefficients 

7fc a coefficient in the recursion of the wave function coefficients 

Sk a coefficient in the recursion of the wave function coefficients 

Ak the proportion relation of Ck to Pk{E) 

AV{x) free function in Hq 

Ek a coefficient in the recursion of the wave function coefficients 

e scaled E 

K Cramer's determinant 

^ the group variable 

■0 the wave function 

■0 the pre-exponential function 

^ the initial value solution of the Schrodinger equation 
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Chapter 1 
Introduction 



In this chapter the basic properties of the QES potentials are described 
and the simplest example (the sextic potential) is analyzed. Quasiex- 
actly solvable (QES) potentials are distinguished by the fact that a 
finite number of wave functions and energy levels can be analytically 
calculated [1]. Hamiltonians corresponding to QES potentials can be 
constructed from bilinear combinations of Lie group generators. Hamil- 
tonians arc second order differential operators, and does not have a 
first order derivative term, whereas a bilinear combination of Lie group 
generators does. By transforming the Hamiltonian into a quasigaugc 
invariant Hamiltonian, the two representations (Hamiltonian and Lie 
group generators) of the Schrodinger equation can be made equiva- 
lent through a change of variables. By choosing the basis of the wave 
functions as the direct sum of the representation of the Lie group and 
an arbitrary set from the orthogonal space, the Hamiltonian becomes 
block- diagonal with at least one finite block. The solutions correspond- 
ing to this block may be calculated exactly. This work mainly deals 
with one dimensional QES potentials which arc generated by the s/(2) 
Lie group generators. This work is within the Lie-algebraic formalism 
[2]. An analytic formalism is summarized in the review article [3]. 
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1.1 The Quasigauge Symmetry 



The Schrodinger equation (defined for example in [4]) can be solved in 
numerous ways, exact and approximate. This equation is gauge invari- 
ant even though the Hamiltonian is not. A gauge invariant Hamiltonian 
can be created by replacing the derivatives by covariant derivatives. 
Quasigauge invariance is defined as invariance under imaginary pliase 
gauge transformations. The quasigauge Hamiltonian is defined as a 
quasigauge invariant operator which is derived from the Hamiltonian 
by a real correction to the derivative (an imaginary gauge field). The 
Schrodinger equation is not quasigauge invariant but if the Hamiltonian 
is replaced by the quasigauge Hamiltonian we can achieve a quasigauge 
invariant formalism. 

For all potentials, the solution of the Schrodinger equation 



HiP{x) 



1 



+ V{x) 



ip{x) — Eip{x) 



can be written as the product 

V'(x) 



(1.1) 



(1.2) 



where a{x) is the imaginary phase, e""*^^-* is the exponential factor 
and ip{x) is the pre-exponential function. Substituting ip{x) in the 
Schrodinger equation leads to a new equation for the pre-exponential 
function ip{x) 



Ho<l> = 



(1.3) 



where Ai{x) = -^a{x). The operator Hq is called the quasigauge 
Hamiltonian and the new equation is invariant under the quasigauge 
transformation 

^^^e-"^-), Ai^A-^. (1.4) 

OXi 

The eigenvahies (energy levels) of equation (1.3) and of the Schrodinger 
equation (1.1) are identical. The eigenfunctions of the two equations 
are related through (1.2). 
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1.2 The Quasigauge Hamiltonian 



QES potentials are distinguished by the fact that their quasigauge 
Hamiltonian may be written as a block diagonal matrix with at least 
one finite block. 





hi,2 


CO 




\ 


h2,i 


h2,2 




hN-2,N 









hN-l,N-l 


hN-l,N 






hN,N-2 


hN,N-l 


hN,N 







hij 7^ / 



(1.5) 



The energy levels and wave functions corresponding to the finite block 
can always be calculated. Quasigauge Hamiltonians can be generated 
by taking bilinear combinations of Lie group generators 

HG^Y.^^abT^'T' + Y.CaT'' (1.6) 

a,b a 

in a basis {\tp)} which is the direct sum of the jth representation Rj of 
the Lie group and an arbitrary set from the orthogonal space 

{IV')} — {elements of Rj + arbitrary set from orthogonal space}. 

(1.7) 

A first order differential representation is chosen for the generators 
so the quasigauge Hamiltonian is a second order differential operator 
equivalent to (1.3) up to a change of variables. 

Many examples of QES potentials are found in the literature: one 
dimensional potentials [2], one dimensional radial potentials [5], QES 
potentials in higher dimensions, spatial and Grassmann (supcrsymmct- 
ric) [6, 7, 8, 9, 10, 11]. Many particle systems of QES potentials are 
analyzed for example in [12, 13]. Because the QES potentials are an 
intermediate class between the exactly solvable potentials and the non 
solvable potentials they were used as a tool in many calculations, for 
example [14, 15]. 
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1.3 One-Dimensional QES Potentials 



This work concentrates on one-dimensional QES potentials. This is be- 
cause in one-dimensional quantum mechanics the ordering of the energy 
levels and wave functions is known (for a review of oscillation theorem 
see appendix A). The one-dimensional QES potentials discussed in this 
work are related to the s/(2) Lie algebra [1] defined by the commutation 
relations 

[r+, T-] = 2r° , r°] = -r+ , [r-, r°] = (i.s) 

This algebra may by represented by 2 x 2 matrices that are non-singular 
under matrix multiplication with determinant [16, 17]. The s/(2) 
group representations are characterized by a single parameter j, where 
the jth representation is of dimension 2j + 1. The differential represen- 
tation of the generators corresponding to j is 

r+ = m-e^^ (1-9) 

= -i + C;^ (1-10) 

r- ^ A. (i.u) 

The jth polynomial representation is composed of powers of the group 
variable, from zero to 2j 

= (1-12) 

The basis of the solution space is ■C^-'j "^2^+21 "^2^+31 •••} where 

{'^2j+k} is an arbitrary set that spans the orthogonal space. The pre- 
exponential functions il){x) of the QES sector are polynomials of degree 
2j and the general wave function, which solves the Schrodinger equation 
(1.1), has the form 

iP{x) = e-"(^)vi = e-"(^) Ck^\^)- (1-13) 

fe=0 
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Writing Hq as an explicit expression in the variable ^ results in 

a,6=±,0 a=±,0 
= -^P4(0^+P3(0|+P2(0 (1.14) 

where Pi{i) is a polynomial of degree i. Without loss of generality 

the coefficients Cah can be taken to be symmetric {cah = Ci,a). The 
polynomials P4(^), PsiO^ P^iO under this assumption are given by 

P4(0 = -2[c++e'-2c+oe'+(coo-2c+_)e' + 2co-e + C— ] 

PsiO = -2(2j-l)c++e' + (3(2i-l)c+o-c+)f 

+ (-(2j - l)(coo - 2c+_) + Co)e + (-(2j - l)co- + c_) 

P2(e) = 2j(2j-l)c++f + (-2j(2j-l)c+o + 2jc+)e (1.15) 

The transformation between the Schrodinger picture and the quasi- 
gauge picture is found by comparing the quasigauge Hamiltonian in ^ 
(1.14) to the quasigauge Hamiltonian in x (1.3). 

Hg--1[-^-Ax)J + V{x). (1.16) 

The change of variables x{^) is found by comparing the coefficients of 
the second derivatives 

Comparing the coefficients of the first derivatives yields the quasigauge 
field A{x) 



/Pi 

and the exponential behavior a{x) 



:i.i8) 
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AV{x) is defined as 

AVix) = p,iOk^^i.) (1-20) 
The free functions are compared to give tlie Sclirodinger potential V{x) 

V{x)^AV{x)+^-A'{x)-^-A'{x). (1.21) 

It is now possible to take any bilinear combination of s/(2) generators 
and find the corresponding QES potential [1, 2]. 

The Schrodinger equation of any QES potential can now be replaced 
using (1.13) to get a set of 2j + 1 equations for the wave function coef- 
ficients {Cfcl^-^Q. This set of homogeneous equations has a non-trivial 
solution if the determinant of the coefficient matrix (the secular poly- 
nomial) is zero. This is the quantization condition of the energy levels. 
The generators can be defined differently [18] so that the Hamiltonian 
is a linear combination of generators but that does not change the range 
of QES potential obtained by this formalism and their solutions. The 
properties of the Bender-Dunne polynomials and other sets of poly- 
nomials that are the wave function coefficients of the QES potentials 
are analyzed in many articles. The orthogonality of the wave function 
coefficients is discussed for example in [19, 20]. 



1.3.1 Example 1 - The Sextic Potential 

The simplest and best known example of a QES potential is the sextic 
potential with a quantized coefficient before the harmonic term. 

V{x) = ^x''-^i8j+3)x\ (1.22) 

This example has been analyzed in many articles. 

This potential has an a; — >■ —x symmetry so a natural choice for a 
new variable is ^{x) = x^ . In ^ the kinetic term of the Hamiltonian is 

~2rf^2 =- 'de'di ^ ' ' 

which can be written in the generators of sl{2) (1.9-1.11) as 

= -2T^T- - (2j + 1)T-. (1.24) 
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Figure 1.1: The sextic potential with J = §• 

The potential energy term can be any linear combination of the genera- 
tors. The sextic potential (1.22) is the result of choosing the quasigauge 
Hamiltonian 

Hg = -2r°r--(2j+i)r--2r+ = -2e^+(2e'-i)^-4je. (1.25) 

Using the definitions of equations (1.17-1.19) the wave functions have 
the form 

2i 

V;(a;) = e-"'/^^Cfex2*= (1.26) 

fc=0 

In [21] the following recursion relation for the coefficients of the wave 
functions 

= 2(2j + 1 - k)Ck-i + ECk + {k + l){2k + l)Ck+i (1.27) 
with the boundary conditions 

C_i = 0, C2j+i = (1.28) 
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Figure 1.2: The wave functions of the sextic potential with j = |. 



is derived. In [22] solutions of the form 

k=o \'^'^)- 

are substituted in the Schrodinger equation to give the recursion rela- 
tion 

Pk+,{E) = EPk{E) - 2k{2k - l)(2i + 1 - k)Pk-i{E). (1.30) 
For A; = 2_7 + 1 the last coefficient of the recursion vanishes and 

{E) = P2,+i{E)Qr,{E) (1.31) 

which means that all the polynomials P2j+i+niE) have the same zeros 
as P2j+i{E). This leads to two conclusions. First, for all the zeros of 
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P2j+i{E) the wave function (1.29) is normalizable because the sum is 
truncated. Second, the polynomial P2j+i{E) is the secular polynomial 
of the set of equations (1.27) if the initial conditions are chosen as 

Po = 1, Pi = E, (1.32) 



And the initial value solution of the Schrodinger equation (1.29) is the 
generating function of the secular polynomials. The first polynomials 



J 








E 





E 





1 

2 


E^- 


2 




±^2 


1 


E'- 


16E 




0,±4 


3 
2 


E^- 






±1.77966, ±7.53875 


2 


E^- 


160^3 + 2944^ 




0, ±4.60568, ±11.7808 


5 
2 


E^- 


350E^ + 2071QE^ - 


81000 


±2.05066, ±8.35389, ±16.6135 


3 




672E^ + 95616^3 - 


2045952E 


0, ±5.09375, ±12.7813, ±21.9703 



Table 1.1: The secular polynomials and QES sector energy levels of the 
sextic potential. 



and energies are hsted in table (1.1). 

The sextic potential has an additional symmetry, known as the en- 
ergy reflection symmetry or self duality. This symmetry was first an- 
alyzed and used in [21] to explain the large j approximation of the 
highest energy level of the QES sector E^ = ^ \/| the reflection of 
the lowest energy at the bottom of the potential well. More potentials 
were found to be self dual. Other potentials were found to satisfy a 
weaker form of the symmetry which is known as duality [23, 24, 25]. 

Parity of the Wave Functions 

The solutions (1.29) are even functions. The odd wave functions of 
the potential (1.1) arc not a part of the quasiexactly solvable sector. 
A very similar potential with an odd QES sector can be derived by 
transforming the quasigauge Hamiltonian (1.25) by: 

Hg^Hg-2T-. (1.33) 
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This new quasigauge Hamiltonian is related through (1.21) to the po- 
tential 

V{x) = ^x'-^i8j+5)x', (1.34) 

which is similar to the original potential (1.1) but has odd QES wave 
functions. The wave functions of this quasigauge Hamiltonian have the 
form 

^(a;) = e-^*/^x ^ Cfex''^ (1.35) 

fe=0 

and the coefficients Ck are connected by the recursion relation 

= 2(2j + 1 - k)Ck-i + ECk + {k + l){2k + l)Ck+i (1.36) 

which is different from (1.27) but has the same initial conditions. The 
even wave functions of this potential are not a part of the quasiexactly 
solvable sector. 



1.4 Outline 

This work concentrates on one dimensional QES potentials. 

In chapter (2) the QES potentials and the various properties of their 
exact solutions are discussed. The wave function coefficients are given 
a new explicit form in terms of the coefficient matrix. A new formula 
for calculating the secular polynomials is derived. Several special cases 
are discussed including supersymmetric potentials. 

In Chapter (3) of this work, a new technique based on the WKB 
approximation is developed, and is used to approximate the highest 
energy levels for large j. This technique combines the initial value 
WKB approximation with the saddle point approximation. It applies 
to all cases that have a generating function of the secular polynomials 
that is the initial value solution of the Schrodinger equation. Two 
explicit examples of using this technique are given. 

In chapter (4) of this work the results and conclusions are summa- 
rized and discussed. 
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Chapter 2 
Exact Solutions 



QES potentials are distinguished by having a finite number of exact 
solutions. The number of the exact solutions is potential dependent. 
This attribute is related to the symmetry of the Hamiltonian. In order 
to explicitly see the underlying symmetry it is more convenient to work 
with an alternative operator: the quasigauge Hamiltonian, Hq, which is 
symmetric under quasigauge transformations. General exact solutions 
are found and the methods are applied to several one-dimensional and 
one-dimensional SUSY potentials. 



2.1 General Solutions 



For the one-dimensional QES potentials it is possible to replace the 
Schrodinger equation with the equation 



a,b=±,0 a=±,0 



(2.1) 



where T", a = ±, are the generators of the s/(2) Lie group. In the jth 
polynomial representation the solutions ip are given by a power series 
in the group variable ^ 

i^ = T.Cke (2.2) 

k 

Equation (2.1) can now be transformed into a recursion relation for the 
coefficients Ck- For the most general one-dimensional case this relation 
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is a 5 term recursion relation: 



.0 = -c++{2j + 2-k){2j + l-k)Ck-2 

- [c+o{2k-2j-l) + c+]{2j + l-k)Ck-i 
+ [E + (coo - 2c+_)(2j -k)k- cok] Ck 

- [co-(2A;-2j + l) + c_](A; + l)Cfe+i 

- c__(A; + l)(A; + 2)Cfe+2. (2.3) 

If we define a^, (3k, 'jk, S^, to be the functions 

ak = -c++{2j + 2-k){2j + l-k) 

Pk = -[c+o(2A;-2j-l) + c+](2j + l-A;) 

7fc = ^+[(coo-2c+_)(2j-/c)-co]/c 

5k = -[co-(2A;-2j + l) + c_](A; + l) 

€k = -c__(A; + l)(A; + 2) (2.4) 

then the recursion relation can be written as 

= akCk-2 + PkCk-i + IkCk + 5kCk+i + (2.5) 

Wave functions arc boundary value solutions of the Schrodinger equa- 
tion. This remains unchanged though it is now a difference equation 
and not a differential equation. In order to truncate the sum (2.2) 
and ensure the normalizability of the Schrodinger wave function the 
boundary conditions are chosen to be 

C_2 = 0, C_i = 0, C2J+1 = 0, C2,+2 = 0. (2.6) 

This, with the special form of the coefficients (2.4) results in the nor- 
malizable wave function 

^{x) = e-''^^^J2Ckeix), (2.7) 

fe=0 

which has no more than 2j + 1 zeros. 

The set of equations giving the coefficients {Ck}^=Q can be written 
in matrix form as 

M2j+iC = (2.8) 
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where 'M.2j+i is the (2j + 1) x (2j + 1) coefficient matrix: 





7o 


5o 


eo 













\ 




Pi 


71 


5i 
















0L2 




72 



















Qi3 






<^2j-3 


e2j-3 

















/32,-2 


72j-2 


<^2j-2 


e2j-2 














q;2j-i 




72J-1 






I 













«2j 


P2j 


l2j 


J 



(2.9) 



(oifc, 7fc, efe as defined in (2.4)), and C is the vector of variables 
C: 

( Co \ 

C2 



c 



(2.10) 



C2i-2 
C2J-I 

V C2, y 

The orthogonahty of the wave function coefficients as functions of the 
energy E is analyzed in [19, 20] and will not be addressed in this work. 
This set of homogeneous equations has a non trivial solution if and only 
if the determinant of the coefficient matrix vanishes 



P2,+i(^) = detM(2,+i) = 0. (2.11) 

This is the quantization condition for the energies of the system. The 
polynomial P2j+i{E) is known as the secular polynomial (or secular 
determinant). 



2.1.1 The Wave Function Coefficients 

The wave function coefficients can be given a general form in terms 
of the minors of M2j+i. This calculation may be generalized to give 
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a solution for a set of n homogeneous equations in n variables. The 
matrix representation of the set of equations has the general form 



M2j+1 C2j+1 



v 



7o 


'^0 


eo 











Co 


A 


7i 


Si 











Ci 




P2 


(32 j -I 


S2j-2 








C2j-2 







72i-l 


S2j-1 






C2J-I 







a2j 


/52, 


l2j 


J 


V 


C2j 



= 0. 



(2.12) 

This set has a non-trivial solution if the determinant of the coefficient 
matrix vanishes (2.11). In such cases the system of equations is lin- 
early dependent and there are less linearly independent equations than 
variables. Since this is a one dimensional quantum system there is no 
degeneracy in the energy levels, so we can expect that there is only 
one free parameter and one extra equation. For a specific QES poten- 
tial the coefficient C2j can never vanish because it will give a solution 
with less zeros than is allowed by oscillation theory (see appendix A). 
This makes C2j a natural choice for the free parameter. The last equa- 
tion is chosen to be redundant. The system of the first 2j equations 
(0 < A; < (2j — 1)) of the system (2.12) for the variables Co, C2j is 
an inhomogeneous system given by 



eo 
Si 



\ 



/ 7o So 

Pi 71 

"2/52 '■• S2j-2 

V ■ • . p2j-l I23-I ) 



( Co \ 

Ci 



/ 





e2i-2C2j 
V S2j-lC2j J 



(2.13) 



V C2,-l J 

which can be written in matrix notation as 

m2jC2,^b2,. (2.14) 
The solution of this inhomogeneous system is given by Cramer's rule: 

dct (Kk) 



Ck 



det ni2j ' 



(2.15) 
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where is the determinant of a matrix that is obtained from ni2j by 
replacing the k^^ column with the vector b2j: 

We have chosen — 1 in (2.15), changing this choice amounts to a 
multiplication of the wave function by a constant and does not change 
the following conclusions since none of these solutions are normal- 
ized. We want to express the solution (2.15) in terms of the minors 
of the original matrix (2.9). It is easy to see that Hk is identical to 
(— l)^-'~'^detM(2j,fc) where M(^n,k) is the minor obtained by eliminating 
the n + 1 row and the k + 1 column of (2.9). The solution for Co, C2j, 
in terms of the minors of M2J-1-1, is 

^ (-l)^^-MetAf,.,,, 

The free parameter C2j can also be included by this determinant ratio 
since 

_ _ detM(2j,2,-) . . 

It can be proved that the boundary condition 6*2^+1 = is satisfied 
automatically since the form of the recursion relation (2.3) for k = 2j — 1 

= a2j-lC2j-3 + P2j-lC2j-2 + {E + 72j-l)C2j_l 

+ S2j-lC2j + €2j-lC2j+l (2-19) 

leads to 

C2j+1 — [(X2j-lC2j-3 + P2j-lC2j-2 

^2j-l 

+ {E + ^2j-i)C2j-i + S2j-iC2j]. (2.20) 
Using (2.17) and (2.18) results in 

(_l)2i+l 

C2j+1 ^ 7— T7 \a2j-lM2j,2j-3 - P2j-lM2j,2j-2 

e2j-i det M^2j,2j) 

+ {E + j2j-l)M2j,2j-l-S2J-lM2J,2j], (2.21) 
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which is equivalent to to a determinant of a matrix with two identical 
rows, automatically giving the boundary condition 



C- 



2j+l 



= 0. 



(2.22) 



It is possible to show that the boundary condition C2j+2 = is equiv- 
alent to the quantization condition (2.11). The recursion relation for 
k — 2j is 

= a2jC2j-2 + /32jC2j-l + {E + 72j)C2j + 52jC2j-l + e2jC2j+2: (2.23) 

and therefore C2j+2 is 

C2j+2 = —[a2jC2J-2 + P2jC2j-l + {E + j2j)C2j]. (2.24) 

e2j 

Substituting the coefficients Ck from equations (2.17), (2.18) and (2.22) 
leads to 



which means that the quantization condition coincides with the bound- 
ary condition C2j+2 = for the most general QES potential. 

These calculations are not valid if det M{2j,2j) = or if = 0, but 
for all other QES potentials the wave function coefficients are given 
by the minors of the coefficient matrix M(2j+i) with the quantization 
condition coinciding with the boundary condition C2j+2 = and the 
boundary condition C2j+\ — automatically satisfied. The special case 
of ejfc = is discussed in section (2.2). 

2.1.2 The Secular Polynomials 

In this section j is kept as a free parameter. The secular polynomials of 
the QES problems which transform into matrices (2.9), which are four- 
diagonal or less, may be calculated by a recursion relation. In these 
cases the recursion for the secular polynomials is of the same degree as 



[a2jM2j2j-2 - P2jM2J,2j-l + {E + ^2j)M2j,2j] 



e2j det M(2j- 2i) 
(-1)^^+^ detM2,+i 
e2j det M^2j,2j) 




(2.25) 
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the recursion for the wave function coefficients (2.23). This does not 
imply that the two sets of polynomials are of the same degree. For 
example for — the determinant of the matrix M^+i as defined in 
equation (2.9) can be calculated by 

Pk+i{E) = {E + ^k)Pk{E)-f3k5k-iPk-i{E) 

+ PkPk-iSk-2Pk-2iE) (2.26) 

with initial conditions Pq = 1, Pi = 70, P2 = 7o7i — Pi^o- Pk{E) is a 
polynomial of degree k whereas the coefficient is not (the case where 
efc = is described in section (2.2)). 

For the most general potentials, which have a five diagonal coeffi- 
cient matrix (2.9), a new recursion relation for calculating the secular 
polynomials is derived. This recursion relation has more than 5 terms 
as shown in the following calculation. 

We start with the coefficient matrix Mfe+i 

/ o/„ A„ c„ n n n ■■ \ 



M 



fc+i 



7o 


4 

















7i 




ei 












^2 


72 

















/93 




4-3 


^k-3 













l3k-2 


7fc-2 


4-2 


^k-2 










Qife-i 


h-1 


7fe-i 


4-1 













Oik 


Pk 


7fe 



(2.27) 



V 

Now let us define the matrix that is a specific minor of M^+i - the 
minor of the term 



/ 


7o 


4 


eo 

















Pi 


71 


4 


ei 














012 


P2 


72 


4 

















as 


/93 




4-4 


efc-4 

















/5fc-3 


7fc-3 


4-3 

















ak-2 


A-2 


7fc-2 


^k-2 




V 













Oik-l 


Pk-l 


4-1 


) 



(2.28) 



21 



The determinant of M^+i in terms of its minors of the forms Mj and 
Ni is 

det(Mfe+i) = (£; + 7fe)det(Mfe) -/3fedet(Nfc) 

+ ak5k-i det(Nfe_i) - ak{E + -fk-i)€k-2 det(Mjfc_2) 
+ Q;jfcQ;fe_iefe_2efe-3det(Mfe_3) (2.29) 

whereas the determinant of Nj^ in terms of its minors of the forms Mj 
and Nj is 

det(Nfc) = 4-idet(Mfc_i) -/3fe_iefc_2det(Mfe_2) 

+ afc_iefc_2det(Nfe_2). (2.30) 

The determinant of the matrix is a polynomial of degree k that is 
denoted hy Qk 

Qk = det(Nfe). (2.31) 
The determinant Pk+i is therefore 

Pk+i = (E + -fk)Pk - PkQk + ak5k-iQk-i 

- 0£k{E + 7fc_i)efe_2Pfe-2 + akak-i€k-2€k-3Pk-3 (2.32) 

and the determinant Qk is 

Qk — ^k-lPk-l ~ /3fc-lCfc-2-Pfc-2 + 0;fe_iefe-2<5fe-2- (2.33) 

To obtain a recursion that will include only the secular polynomials 
the polynomials Qk,Qk-i,Qk-2 must be replaced. A third equation, 
equation (2.32) with A; —> A; — 1 is added 

Pk = {E + -fk-i)Pk-i- /3k-iQk-i + ak-iSk^2Qk-2 

- ak-i{E + 7fc-2)efc-3-Pfc-3 + ak-iak-2ek-3ek-4Pk-4- (2.34) 

The equations (2.32, 2.33, 2.34) form a set of three equations in the 
variables: Qk:Qk-i,Qk-2- Solving these equations and substituting 
Qk-i,Qk-2 in equation (2.34) gives the wanted recursion relation for 
Pk{E) 

= [—ak-2Sk-4Sk-3 + (^k-3(^k-2^k-A\Pk 
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+ [Pk^3{-Pk-2{E + 7fe-i) + ak-iSk-2)ek-4 

+ ak-2Sk-A{{E + 7fc^i)(5fc_3 - Pk-i^k-4.)]Pk-i 

+ [6k-2if3k-3iPk-2f3k-i - ak-i{E + 7A-2))efe-4 

+ Q;fc_25fc-3(-/5fc-l5fc-4 + ak-iek-4)) 

+ Oik^2f3k-l{{E + 7fe-2)4-4 - /3fc-2efc-4)efc-3]-Pfc-2 

+ [afc_2afc_i(5fc_4(5fc_3((5fe_3(5fe_2 - {E + 7fc_3)efc_4 - (-E + 7fc_2)efc_3) 

+ Pk-2{0ik-2Pk-l{E + 7fc_3) + /3fc-3(-/3fe-2/3fe-l 

+ afe-i(£^ + 7fc-2)))efc-4efc-3]^fe-3 

+ [«fc-2ttfc-l/3fc-3(-'5fc-45fc-34-2 + (-E + 7fc-3)5fc-2efc-4 

— Pk-2^k-4^k-3)^k-4 + Oik-2Sk-4{Pk~3Pk-2Pk-l 

+ ak-2{-Pk-i{E + 7fc_3) + afc-i4-3))efc-4efc-3]A-4 

+ [ak-30^k-2{o^k-i{E + 7fc-4)4-34-2 - afc-i/3fc_34-2efc-4 

+ Pk-i{-f3k-2{E + 7fc-4) + afc-24-4)efc-3)efc-5efc-4]-Pfc-5 

+ [Q;fc_4Q;fc_3afc_2(-Q;fc-i4-34-2 

+ (3k-2/3k -l^k-3)^k-ei^k-5^k-4 ]Pk-6. (2.35) 



Initial conditions for this set of equations are the first six polynomials 
Po, Pi: -P2, -P3, -P4, -Pg- This is a 7 term recursion but it becomes a four 
(or less) term relation if one or more of the coefficients 0;^, 5k, 
are chosen to be identically zero. 



2.1.3 Parity of the QES Sector Wave Functions 

For every even potential (not only QES) the Hamiltonian can be sep- 
arated into two blocks with different parities. In the QES potentials 
either one or both of these blocks are QES . If ^(,t) is an odd function 
both blocks are always QES . On the other hand if ^(x) is an even 
function the generic case has only one QES block. In a certain family 
of QES potentials with even ^{x) for some of the potentials the QES 
block is the even one, and for others it is the odd one. The difference 
is in the parity of e"^^^ . 
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QES Sectors of different peirities 

The general transformation from a potential with even parity wave 
functions to a potential of the same family (the same ^{x)) with odd 
parity wave functions can be derived by changing the linear terms in 
the quasi-gauge Hamiltonian (1.14). A family of potentials with even 
^{x) and even e~"^^^ is considered. What is the closest family with the 
same C{x) but odd e-«(^) ? The new Hamiltonian is connected to the 
original by 

H^^'' ^ Hg+ CaT" (2.36) 

a=±,0 

= -\p^i0^2 + b3(0 +Psm I + \P2iO +P2] , (2.37) 

where psiC) and P2{C) are the corrections to psiC) and P2{0 

P3(0 = -c+e + Co^ + £- (2.38) 
MO = 2jc_e. (2.39) 

The odd wave function is 

^om{x) = e-^o'^^^^) J2 Ck^ix). (2.40) 

fe=0 

The coefficients Ck differ from the wave function coefficients Ck in the 
original potential but they are calculated using the same formalism. 
The new relevant functions (1.18-1.20) are 



'Pi ^/Pa 
I P'j^+Ps+Ps ^^ ^ ^^^^ ^ r^(^) P3 ^2.42) 

J P4 J Pi 



(2.41) 



aoddyx) 

^Vom{x) = P2(0+P2(0 = Ay(x)+p2(0l^=^(.) (2.43) 
The wave functions of the QES sector have the explicit form 

Pi 



^Ux) = exp |- /^^"^^ ^dA e-<^^ ^ke{x). (2.44) 
{ J Pi \ t^o 
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The function exp 

{- f^""^ ^d^} is chosen to be odd. As in (1.21) the 



corresponding potential is 



Vom(x) = AVoUx) + IaIm(x)-Ia',Ux) (2.45) 




= nx)+\p2 + A^-b'3 + ~(p^ + lp'4)\ (2.46) 
[ VpI 2 2p4 V 2 /J^^^(,) 



The potential Voddi^) is the most similar to V{x) that has an odd QES 
sector. 

2.1.4 Example 2 - Potential Built From Hyper- 
bolic Functions - Sine 

The hyperbolic sine potential is an example that demonstrates many 
of the things discussed earlier. The quasigauge Hamiltonian in this 
example is 




(2.47) 



(2.48) 



and the relevant functions (1.17-1.20) are 

^{x) = cosh,x 
A{x) = a sinh x 
a{x) = a cosh x 



(2.49) 
(2.50) 
(2.51) 
(2.52) 



AV{x) = —2jacosh.x. 
Equation (1.21) gives the hyperbolic potential 




(2.53) 



The wave function is 



ijj — e 



—a cosh a; 



^Cjk(cosha;)^ 



(2.54) 
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Figure 2.1: The hyperbolic sine potential with a = 1, j = |. 

The recursion relation that corresponds to the Schrodinger equation is 

1 







= a{2j + l-k)Ck-,+ [^E + -k'jCk 



(2.55) 



with the boundary conditions C_i = 0, C2j+i = 0, C2j+2 = 0. The 
recursion relation for the secular polynomials is 



k+l 



Pk - a^k{2j + 1 - k)Pk-i 



-a^k{k - l)(2j + 1 - A;)(2j + 2 - A;)P,_2, (2.56) 



with the initial conditions Pq = 1, Pi = E, P2 ^ + - 2j d/- . 

The odd potential that is closest to this potential is obtained by 
transforming 

Hg^Hg- T°. (2.57) 
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Under this transformation ^(x) is unchanged and the relevant functions 
(2.41-2.43) are 



Aodd{x) — a sinh x — coth x 
'^ (•^)odd — a cosh x — ln(sinh x) 
^V{x)odd = -2jacoshx. 



The odd parity QES sector potential (2.44) is 

Vodd{x) = y sinh^ x-a (2j + coshx + ^. 



The odd wave function is 

V' = sinhxe-""°'''^ 

The recursion relation for the coefficients is 




2i 

Y.Ckicoshx)'' 

k=0 



(2.58) 
(2.59) 
(2.60) 



(2.61) 



(2.62) 



= a{2j + 1 - k)Ck-i +^E + -k{k + 2) 
+ a{k+l)Ck+i-^{k + l){k + 2)Ck+2: 



with the boundary values C_i = 0, C2j+i = 0, 6*2^+2 
recursion relation for the secular polynomials is 



(2.63) 
0. The 



fe+i — 



E+-k{k + 2) 



Pk - a^k{2] + 1 - k)Pk-i 



1 



a\{k - l){2j + 1 - A;)(2j + 2 - k)Pk-2 (2.64) 



with the initial conditions Po = 1, Pi = E, P2 = E"^ + E - 2j a^. 
The wave function is the same as in (2.54) but the coefficients are 
different. 



2.2 Special Case - = 

In all the cases with = (e^ is defined in equation (2.4), ejt = is 
equivalent to c = in Hq) , can be polynomials of degree k if the 
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free parameter is chosen to be Cq = 1. In this section it is proved that 
Pk and Cfc are identical up to an overall proportion factor A^. This 
means that the initial boundary solution of the Schrodingcr equation 
(for a general discussion see Appendix B) is the generating function of 
the secular polynomials. This property was shown in [22, 19] to hold 
for specific examples, and is here expanded to include the general group 
of potentials with = 0. 

If = 0, then the recursion relation for the wave function coeffi- 
cients is 

= akCk-2 + PkCk-i + {E + -ik)Ck + kCk+i. (2.65) 

Since this is a four term recursion relation, only 3 boundary value 
conditions are necessary 

C_2 = 0, C_i = 0, C2j+i = 0. (2.66) 

The recursion relation for the secular polynomials, derived indepen- 
dently from the matrix (2.9) using the known properties of determi- 
nants is 

Pk+i{E) = {E + -fk)Pk{E)-(3kSk-iPk-i{E) 

+ ak6k-i6k-2Pk-2, (2.67) 

with the initial value conditions 

Po = 1, Pi= 70, P2 = 7o7i - A^o- (2.68) 

From the two recursion relations a connection between the wave func- 
tion coefficients and the polynomials can be derived: 

Ck = (-1)' (n Pk{E) = A,P,{E) (2.69) 

The quantization condition P2j+i{E) = coincides with the boundary 
condition C2J+1 = 0. The wave function can now be written as 

M^) = e""^'^ E A,P,{E^)e{x) (2.70) 

k=0 
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where En are the eigenvalues of the Schrodinger equation (2.69) into 
equation (2.23) but leaving the value of E undetermined gives a gener- 
ating function for the secular polynomials of this system. The function 

oo 

*(x) = e-«(-) Yl ^kPkiE)e{x) (2.71) 

fe=0 

where is a free parameter, is the initial value solution of the Schrodinger 
equation (1.1). It is also the generating function of the secular polyno- 
mials Pk{E). 



2.2.1 Example 3 - Potential Built From Hyper- 
bolic Functions - Cosecant 

If the polynomial P4,{^) (1-15) has a double root in ^{x = 0) then the 
function ^{x) is defined separately for the different sections of the real 
line. The Hamiltonian, the potential and the wave functions are chosen 
to have a well defined parity, either even or odd. 

V-{-x) = ±V+{x), x>Q. (2.72) 

The hyperbolic cosecant potential demonstrates this special quahty dis- 
cussed in [26]. It is generated by the quasi-gauge Hamiltonian 

Hg ^ -- {T+T+ + r°T°) + 2 j (T+ + T" ) - jT° 

= eie + 1)^ + [{2j - i)e - ^je - le + 2j]| 

- j(2j - l)f + (2.73) 
The relevant functions (1.17-1.20) are 

e(x) = 5— = -cschx (2.74) 

^ ^ smhx ^ ' 

= _ 2i(csch3x + csch^-l) 

cothx cschx 

a(x) = 2j(— 2a; + sinhx + Incoshx — Insinhx) (2.76) 

^V{x) = -j{2j-l)e + 4j% (2.77) 



29 



and the potential (1-21) is 

V{x) — — -sech^x(4 — llj + 12j cosh 2x — j cosh Ax 

+ (9 + 32j)sinha; + sinh3a;). (2.78) 
In accordance with [26] the potential and wave functions are defined 




Figure 2.2: The hyperbolic cosecant potential with j — \- 

to be even by replacing x with \x\. The recursion relation connecting 
the wave function coefficients is 

= i(2j + l-A;)(2j + 2-fc)Cfc_2-2j(2j + l-A;)Cfc_i 

+ [£; + y]Cfe-2j(A;+l)Cfe+i (2.79) 

with the boundary conditions C_2 = 0, C_i = 0, C2j+i — 0. The 
secular polynomials are also connected by a four term recursion given 

by 

Pk+i{E) = [E+^-]PkiE)-4fk{2j + l-k)Pk-i{E) 

+ 2f{2j + 1 - k){2j + 2 - k)k{k - l)Pfe_2(^),(2.80) 
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with the initial conditions Pn 



-0 = 1. Pi= E, P2 = E^ + \E + 8f. The 
first polynomials and their zeros are given in table (2.1). The relation 



j 




E 





E 





1 

2 


E' + \E~l 


-1.28,0.781 


1 


E'-i + _ 15^ _ 8 


-5.12,-0.5,3.12 


3 
2 


£;4 + jE'^ - '^E'' - ^E + 729 


-11.48,-4.5,2.03,6.95 


2 


E^ + 15£;4 - ^^E-^ - ^E-' 


-20.36,-11.02,-2.06, 




+13412^; + 35072 


6.20, 12.24 



Table 2.1: The QES sector energy levels of the hyperbolic cosecant 
potential. 



between the coefficients and the secular polynomials is 

1 



C'fc 



The wave functions have the form 



P^E). 



(2.81) 



^„ = e-2^(^^°'^^-2^)(tanhx)2^' V ^^Pfc(£;„)(-cschx)'= (2.82) 



k=0 



where E^ is the eigenvalue of the Schrodinger equation. The 4 analytical 
wave functions of the potential (2.78) with j — 3/2 are illustrated in 
figure (2.3). The generating function of the secular polynomials is given 
by (2.71) 

OO 1 

^(x) = e-'^'(^'"^^-'^)(tanha;)2^' V , , Pfc(P)(-cscha;)^ (2.83) 
where j is a free parameter of the potential. 



2.2.2 Example 4 - Polynomial Potential 

The potential discussed in the introduction is one of a family of poly- 
nomial potentials that have the quasigauge Hamiltonian 

Hg = -2T^T- - {2j + 1)T- - uT+ + piT^ (2.84) 
= -2e;^ + K' + /^e-l);|-2jX, (2.85) 
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Figure 2.3: The wave functions of the hyperbohc cosecant potential 
with j = |. 



and the relevant functions (1.17-1.20) 



A{x) 

a{x) 
AV{x) 



The potential (1.21) is 



X 

-(ra^ + /ix) 

4 2 

UX ^ iix 



8 

-2jvx^. 



8 



X 



(2.86) 
(2.87) 

(2.88) 
(2.89) 



(2.90) 
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(where the scxtic potential is a polynomial potential with u = 2, fi = 
0) . The general recursion relation for the coefficients of the wave func- 
tion is 

= u{2j + 1 - k)Ck-i + {E- fxk)Ck + {2k + l){k + l)Ck+i (2.91) 

with boundary value conditions C_i = 0, C2J+1 = 0. The secular 
polynomials satisfy the recursion relation 

Pfe+i ^{E- iik)Pk - iy{2j + l-k){2k- l)kPk-i (2.92) 

with initial conditions Pq{E) = 1, Pi{E) = E. The generating function 
of the secular polynomials is 

*(x) = « ^ O ^ Wr^^(^)^''- (2.93) 

fe=o y-^'^)- 

2.2.3 Example 5 - Exponential Potential 

The exponential potential has no well defined parity because ^(x) has 
no zeros [26] . The quasigauge Hamiltonian for generating this potential 
is 

Hg = -^T^T^-aT+ + cT- + (b-j + ^^T^ (2.94) 
= ~edl + {ae + bi + c)d^-2ajC, (2.95) 



and the relevant functions (1.17-1.20) are 



e(x) = (2.96) 
A{x) = ae" +{h+]^+ ce"^ (2.97) 



2. 

a{x) = ae" +{h+]^x- ce"^ (2.98) 
AV{x) = -2aje'. (2.99) 



The potential (1.21) is 



V{x) = ^[ae- + {h- 2j)f + ^[ce"^ + (6 + 1)]^. (2.100) 
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The wave function coefficients recursion relation is 



= a{2j + l-k)Ck-i- 



Ck-c{k+l)Ck+i (2.101) 



with the boundary conditions C_i = 0, C2j+i = 0. The recursion 
relation for the secular polynomials is 



Pk{E) 



+ ack{2j + l-k)Pk-i{E) 



(2.102) 



with the initial value conditions Pq{E) — 1, Pi{E) — E. The relation 
between and Pk{E) is 



k\c 



(2.103) 



2.3 Special Case - Exactly Solvable Po- 
tentials 

The exactly solvable potentials can be described as quasiexactly solv- 
able potentials with an infinite quasiexactly solvable sector [1, 27]. This 
requires that the quasigauge Hamiltonian Hq docs not depend on the 
dimension of the representation j. The general quasigauge Hamiltonian 
that is independent of j is 

/fG = -^Q2(0^ + Qi(0^, (2.104) 

where (^2(0 is a polynomial of the second degree or less 

Q2(0 = -2cooe' - 4co_e - 2c__, (2.105) 
and (5i(0 is ^ polynomial of the first degree or less 

QiiO = (Co - (2j - l)coo)e + (c- - (2j - l)co-). (2.106) 
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ak and (3^ in (2.4) are zero for all exactly solvable potentials. The 
coefficient matrix (2.9) is therefore top triangular, which means that 
its determinant is equal to the infinite polynomial 

oo 

PiE) = n 7fc- (2.107) 

k=Q 

All of the wave functions can be analytically calculated and are of the 
form 

n-l 

^n^e-''^'^J2Ck{En)e{x). (2.108) 

fc=0 

Unlike the QES potentials, the nth wave function is a polynomial of 
degree n. The coefficients Ck are given by the recursion relation 

0^{E + 72,)C2,- + S2jC2j-i + e2jC2j+2 (2.109) 

For Cfc = or 5fc = all of the coefficients are known in closed form. 
The best known examples of exactly solvable potentials are given in 
table (2.2). More examples can be created by the factorization method 
known as Infeld-HuU [28, 29]. 



V{x) 


Hg 


Harmonic Oscillator 


-\T-T- +ujT^ 


Harmonic Oscillator — \ 


_2r"r- - {2j + l)T- + 2T" 


Morse Potential A{e-'^°^ - 26""^) 


-°^T^T^ + V2AaT- 


Poschl-Teller Potential ^ 

cosh ax 


+(yi + 8^-(2j + i))ro] 



Table 2.2: The exactly solvable potentials and their sl{2) representa- 
tions. 



2.3.1 Example 6 - Harmonic Oscillator 

The Harmonic oscillator can be described by two different quasigauge 
Hamiltonians. The first, — iT~T~-|-a;T° , gives both odd and even wave 
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functions and their energies. The ahernative representation separates 
the problem into even and odd parts but has the initial value solution 
that is used in chapter (3). The quasigauge Hamiltonian of the even 
states of the harmonic oscillator is given by 



Hg = -2T"T- - (2j + 1)T- + 2T" 



with the relevant functions (1.17-1.19) 

X 



A{x) = = 



AV{x) = 0. (2.111) 
The potential (1.21) is 

Vix) = \x' - \ (2.112) 

The quasigauge Hamiltonian (2.110) gives the even wave functions. The 
odd wave functions are given by a correction of — 2T~ to (2.110) that 
has odd e~"^^) but the same potential V{x). The even wave functions 
are ^ 

^e^en ^ ^-Ix^ ^ Cfe(£;^^"")x''= (2.113) 

A:=0 

The wave functions coefficients are given by a two-term recursion 
relation 

= (E - 2k)Cu + (2A; + \){k + l)Cfc+i (2.114) 

with the boundary condition C_i = 0. The closed form of the coeffi- 
cients Cjfc can be calculated 

^^(^) = 7Wn(^-2^)- (2.115) 

l^'^J- n=0 
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The even energy levels that correspond to these solutions are 

^even ^ 2n. (2.116) 

The wave functions for the even states have the closed form 

2i {_o 



The polynomials P2j+i{E) are 



■fc-i 

2k 



n {E^r - 2^) 

.n=0 



(2.117) 



PI]Ti{E)^J{{E-2k), (2.118) 

fc=0 

and the even wave functions are then written as 

^e.e„ ^ ^-ix^ 1^ ^Pr"(£;2,)x^'. (2.119) 

fc=o y^'^y- 

The initial value solution of the Schrodinger equation is the generating 
function of the polynomials P2j^^{E) 

OO / r)\fc 

^even ^ ^-ix^ ^ ^_^pe^en^^^^2fe_ ^2.120) 

fc=o (2^)' 

Now V{x) and £■ can be redefined to be 

V{x)^^x', £;r" = 2n+^, (2.121) 
The form of ^^'^^^ is unchanged. 



2.4 Special Case - A Complete Separation 
of Even and Odd States 

A system which can be separated into even and odd parts is obtained 
from (2.12) by taking = 0, 6k — 0. The recursion is a 3 term relation 
of the form: 

= akCk-2+{E + ^k)Ck + ekCk+2 (2.122) 
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with boundary values 

C_2 = 0, C_i = 0, C2,+i = 0, C2,+2 = 0. (2.123) 

Prom a practical point of view, since ^{x) is an odd function the total 
wave function can be written as a sum of an even function and an odd 
function with independent recursion relations 

^ ^^even ^^odd_ (2.124) 

The secular polynomial of the whole system is a product of the quan- 
tization conditions of the separate systems 

2j [j] b'-^l 

P2j+i{E) = n(^ - = n(^ - m^- Ef") 

i=0 i=0 i=0 

= g[,+i](E)g^.+ i](E). (2.125) 
The wave function is 
V'2,+i(x) = ^DTii^) + KA^) (2.126) 

- p-aW V 9]l^p'^kf , -a(x) 'V-'^ <ik{E) ^2fc+l/ n 

" V (2^^)! ^ ^ V (2A; + 1)!^ ^ ' 

Even States 

Let us define a new integer variable m = 0, 1, [j]; k — 2m. The even 
wave function coefficients Dj^ — C2m satisfy the recursion relations 

= a2mDm-l + {E + -f2m)Dm + e2m^m+l, (2-127) 

with the boundary conditions D_i = 0, D2j+i = 0. The even energy 
polynomials satisfy the recursion relation 

Ql^E) = {E + l2m)QUE) - a2rr,e2m-2QL-l{E) (2.128) 

with the initial conditions Qo = 1, Qi = E. This problem is solved 
by the same rules described in section (2.2). The relation between the 
coefficients and the polynomials is again a simple proportion relation 

m—l 1 

Dm = {-irQmiE) n — (2-129) 

i=0 
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The quantization condition for the even states is 



(2.130) 



Odd States 

The integer m is redefined to be m = 0, k = 2m+l. The 

odd coefficients dm = C2m+i satisfy the recursion relations 

= Oi2m+ldm-l + {E + '^2m+l)dm + ^2m+ldm+l, (2.131) 

with the boundary conditions d^i = 0, d2j+i — 0. The odd energy 
polynomials q„i satisfy the recursion relation 

Qm+liE) = {E + 'y2m+l)qm{E) - Q;2m+ie2m-l?m-l (^) (2.132) 

with the initial conditions qo = 1, qi = E + 'ji. The relation between 
the coefficients and the polynomials is once more a simple proportion 
relation 



m—l 



dm = {-ir n 



1 



-qm{E). 



i=0 ^2i+l 

The quantization condition for the odd states is 



(2.133) 



(2.134) 



2.4.1 Example 7 - Potential Built From Hyper- 
bolic Functions - Cosine 

The hyperbolic cosine potential is an example of complete separation 
of even and odd states. It is generated by the quasigauge Hamiltonian 



(2.135) 



-o(l-r)X+ < - (2j - l)e(l - r) d^-3{2j-l)i\ 



with the relevant functions (1.17-1.20) 
^{x) — tanhx 

A{x) — (a cosh^ X — 2j) tanh x 
-- — cosh^ X — 2j ln(cosh x) 



a[x) 
AV{x) 



—j{2j — 1) tanh^ x. 



(2.136) 
(2.137) 

(2.138) 
(2.139) 
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The potential (1.21) is 

V{x) = ^ cosh^ X - I (a + 4j + 2) cosh^ x. (2. 140) 

The recursion relation for the coefficients has the form of (2.122) 




Figure 2.4: The hyperbolic cosine potential with a = 1, j 



= ^{2j + 2-k){2j + l-k)Ck-2 

+ [E + k{2j-k-a)]Ck + ^(k + l)(k + 2)Ck+2. (2.141) 
The boundary value conditions are given by C_2 = 0, C_i = 0, C2j+i — 

0, C2J+2 = 0. 

The recursion relation (2.127) for the even state coefficients is 

= (j + l-m)(2j + l-2m)D^_i 
+ [E + 2m{2j -2m- a)]Dm 
+ (2m+l)(m+ 1)^^+1, (2.142) 
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with boundary conditions D_i = 0, -Dy+i] — 0. The recursion for the 
even energy polynomials Qm is 

Qm+i = [E + 2m(2j - a- 2m)]Q„, 

- (j + 1 - m)(2j + 1 - 2m)(2m - l)mg^_i, (2.143) 

with initial conditions Qq = 1, Qi — E. The relation between the 
coefficients and the polynomials is 

f— 2)™ 

Dm = ^-^Qm{E) (2.144) 
The even wave function is 

bl n, ( Tpeven\ 

^e.en ^ ^-a(.) ^ ^k[J^n } ^2k^^y (2.145) 
l^'^J' 

The recursion for the odd state coefficients is 

= (2j + l-2m)(j -m)d^_i 

+ [£;+(2m+l)(2j -2m-a-l)]d^ 

+ (2m + 3)(m + l)d^+i, (2.146) 

with boundary conditions (Iq — 1, = 0. The recursion for the 

odd energy polynomials is 

qm+i = [E + {2m + l){2j -2m-a-l)\qm 

- (2i-2m + l)(i-m)(2m+l)mg^_i, (2.147) 

with initial conditions q^ — 1, qi — E — 2j . The relation between the 
polynomials and the coefficients is 

(—2)'^ 

= (^^«-<^> <^-"**' 

The odd wave function is 

^oM ^ ^-ai.) ^^e'^\x) (2.149) 

The odd and even energy levels are interleaved, with the ground 
state always being even. The first polynomials and energy levels are 
given in table (2.3). The lowest energy levels are very close because of 
tunneling between the two minima of the potential. 
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J 




E+ 




E_ 





1 




1 




1 

2 


E 





E — 1 


1 


1 


E^-2E-1 


-0.414214 
2.41421 


E 





3 

z 


E^-3 


-1.73205 
1.73205 


E'^-2E-6 


-1.64575 
3.64575 


2 


E^ - 2E 
-20E + 24 


-4.17226 
1.14399 
5.02827 


E'^ + 2E-9 


-4.16228 
2.16228 


5 
2 


E^ + AE^ - 28E 


-7.65685 


3.65685 


E^ + E^ 
-49E + 15 


-7.65606 
0.308667 
6.3474 


3 


E^ + 4^3 
-102^2 + 12E 
+585 


-12.150119 
-2.294220 
2.715272 
7.729066 


E^ + lOE"^ 
-36E - 120 


-12.150071 
-2.246438 
4.396509 


7 
2 


E^ + 16E^ 
-62^2 _ 528^ 

+945 


-17.645885 
-5.732662 
1.623139 
5.755408 


E^ + 12^3 
-144^2 - 664E 
+2100 


-17.645883 
-5.728016 
2.286225 
9.087674 



Table 2.3: The QES sector energy levels of the hyperbolic cosine po- 
tential. 



2.5 Super symmetric Potentials 

This section deals with supersymmetric (SUSY) potentials that can be 
represented by 2 x 2 matrices or by Grassmann variables. The wave 
functions are described either as two component spinors or as functions 
of two variables, one of which is Grassmann. The super algebra dis- 
cussed in this section is s-u(2/l), which is the simplest super algebra to 
give QES potentials. Other groups with a more complex structure can 
be chosen to give SUSY potentials that have more spinor components 
or more variables [30]. The su{2/l) algebra includes four even gener- 
ators and four odd generators. The even generators T", a = ±, are 
the generators of the su{2) algebra and J is the u{l) generator. The 
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odd generators are and with q = 1,2. The commutation and 
anti-commutation relations of these generators are given by 

[T+,T-] = 2T°, [T^,T°] = ^T^, [J.T"] = 

[T+, Qi] = Q2, [T+, Q2] = 0, [T-, Q,] = Q,, [T~, Q,] = 

[T-, Qi] = Q2, [T-, Q2] = 0, [r+, Q2] = Qi, [T+, Qi] = 

{Qi,Q2} = -T+, {Q2,Qi} = -T- 

^{Qi, Qi - Q2, Q2} = T°, i{gi, Q, + Q2, Q2} = J (2.150) 



These generators can be realized by the operators 

Q 



= -j + ^d^ + ^Ode ^ \Q2 \id^ 



'e 



J ^ -3- \ede [Q2 ) [ -ed^ J 

(2.151) 

^ is a Grassmann variable reahzed by 



T) 



1 




0' = U)=(?)' (2-152) 
or an operator realized by 2 x 2 matrices 

= (7+ = (CTI +i(72)/2 

de = (7- = ((71 - i(72)/2, (2.153) 
where (7j are the Pauli matrices 

-(?j)--(roO'-(J-0- 

Through the same formalism as in the non SUSY case, the Schrodinger 
equation is replaced by the equation 

Hoi) = Ei> (2.155) 
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where Hq is a a 2 X 2 matrix which is a bihnear combination of all 8 
generators denoted here by 

Hg^YI ('abT'^T' + ^ c„T", (2.156) 

a,b a 

and ^0 is a 2-component wave function 

This equation is invariant under the quasigauge transformation (in the 
SUSY potentials it is a matrix and not a scalar transformation). 



2.5.1 Example 8 - The Simplest Supersymmetric 
Potential 

The simplest known example of one-dimensional supersymmetric QES 
potentials is given by the quasigauge Hamiltonian 

Hg = -{T°,T-}-2jT--Q2Qi + aT° 

- 0{Q2T- + Qi + 2jQ^) + '-aPQ2-'-PQi. (2.158) 

As in example (1.3.1), the change of variables is ^{x) = x^. The poten- 
tial is the 2x2 matrix 



V{x) = l(a^-/3V + ^2(-2i/3-^ + ^ 



a I3x^\ I3x^ 

— tan cos 

4 2 / 2 



f ^ (3 a(5x^ a (3x^\ (3x^ 
+ <73 -2j/3-- + ^--cot— sm — 



(2.159) 



The new methods described in section (2.1) are applied to this potential. 
The Schrodinger equation gives rise to wave functions with coefficients 
that satisfy the recursion relations 



= (2A; + 1)(A;+ 1)4+1 



E-a k 



1 



dk 
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= {2k+l){k + l)uk+i + 
+ i(5{2j-k)dk, 



ia(3 



-Uk-l 



E-aik + 



(2.160) 



with the boundary conditions d-i = , U-i = , U2j = , d2j+i — 0. 
It is convenient to define the functions Ak, Bk, Ck, D^, F^, Gk, 
as 



A, — _ 

Bk = E-a(k-{ 



i/3 (2j + l + k 



Gk 
Hk 



i(5{2j-k) 
E-a(k + \ 



(2.161) 



(2A; + 1)(A; + 1). 



Dk = (2A; + 1)(A; + 1) 
The relations (2.160) then can be written as 

= AkUk-i + Bkdk + CkUk + Dkdk+i 
= Fkdk + GkUk + HkUk+i 



(2.162) 



The recursion relations for Uk and dk can be separated to give a 4 term 
recursion for the upper wave function coefficients 

= [—Ak+iFk+iFk+2]uk 

+ [{Bk+iGk+i — Ck+iFk+i)Fk+2]uk+i 

+ [Dk+iFk+iGk+2 + Bk+iFk+2Hk+i]uk+2 

+ [Dk+iFk+iHk+2]uk+3, (2.163) 

and a 4 term recursion for the lower wave function coefficients 

= [Ak+iFk{Ck+2Gk+i — Ak+2Hk+i)]dk 

+ [Ck+2Fk+i{Ck+iGk — Ak+iHk) 

+ Bk+iGk{Ak+2Hk+\ — Cfc+2G'fc+i)]dfc+i 

+ [Ak+iBk+2HkHk+i + Gk{—Ck+2Dk+iGk+\ 

+ Hk+i{Ak+2Dk+i — Bk+2Ck+i))]dk+2 

+ [Dk+2Hk+i{Ak+iHk- Ck+iGk)\dk+z. (2.164) 
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The system of 4j + 1 equations described in (2.160) can be written in 
matrix form. Choosing the ordering of the variables to be 



— 



f do \ 

di 

C?2j-1 
U2j-1 
d2j 



(2.165) 



gives the (4j + 1) x (4j + 1) matrix Mij^i the form 



4i+l 



Bo 
Fo 











Co Do 
Go 

Fi 





Ho 












V 











A 








2i-i 









H2,-2 
B2j-1 C2j-1 D2j-l 





H2j-2 







G 



A 



2j 





^2, y 



(2.166) 



where B, C, D, F, G, are as defined in (2.161). The Schrodinger 
equation can be written as 



M^j+iC = 



(2.167) 



The determinant of the coefficient matrix P^j^i^E) = detM4j+i is the 
secular polynomial of the system and the zeros of this polynomial are 
the energy levels of the potential. A recursion relation connecting the 
secular polynomials Pk{E) can be found using the technique used in 
section (2.1.2). This recursion is the four term relation 

= [Ak-lDk-2 — Bk-lGk-2]P2k+l 

+ [~^fe-i-Bfe-Dfe-2Gfe_i + Cfc_2(^fe£>fe-iFfe-i 
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+ -Bfc(— Cfc-i-Ffc-1 + Bk_iGk-i))]P2k-i 

+ [Ak-i{—AkDk-2Dk-iGk-2 + Bk{Bk-iCk-2Hk-2 

+ Dk-2{Ck-iGu-2- Ak-iHk-2)))Fu-i]P2k-:^ (2.168) 

+ [A.k-2A.k-iDk-2{BkCk-i — AkDk-i)Fk-iHk-3]P2k-5 

with 3 initial value conditions. 

2.6 Summary of Results 

In this chapter the exact solutions of some QES potentials and some 
exactly solvable potentials are analyzed. Solving the most general QES 
potential in one dimension gives a set of recursion relations for the wave 
function coefficients. The wave function coefficients are found to be pro- 
portional to the minors of the coefficient matrix as a result of applying 
Cramer's rule for the calculation of solutions of anharmonic systems of 
equations. This result explains the connection between the boundary 
conditions and the quantization condition, which is the vanishing of the 
secular polynomial. The recursion relation of the secular polynomials 
in the most general case is found to have more terms than the recur- 
sion relation of the wave function coefficients. The different parities 
of QES potentials are compared and the transition from even to odd 
parity potentials defined. For all cases with = it is proved that the 
generating function of the secular polynomials solves the Schrodinger 
equation with initial conditions. The one-dimensional exactly solvable 
potentials are analyzed using the QES formalism and the difference 
between exact solvability and quasiexact solvability is discussed. For 
supersymmetric potentials the recursion relations for the upper and 
lower coefficients are separated, and a recursion relation is found for 
the secular polynomials. 
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Chapter 3 
Approximations 



Typically, the QES potentials have a free parameter j which determines 
the size of the QES block. For j ^ 1 the calculation of the energy lev- 
els is complicated. A method of approximation is needed, which can 
simplify the calculation of the secular polynomials and energy levels 
for large j. In section (2.2) the generating function of the secular poly- 
nomials has been defined. In this chapter the generating function is 
approximated using the WKB approximation with initial value condi- 
tions (and not with boundary value conditions), and using the saddle 
point approximation. A new approach to approximating the energy 
levels is described in this chapter. 

3.1 WKB Approximation of the Generat- 
ing Function 

In section (2.2) it is proved that for all cases with €k — (2.4) the initial 
value solution of the Schrodinger equation is the generating function of 
the secular polynomials (2.71) 



oo 




(3.1) 



k=0 



with initial conditions 



Po{E) = l, P,{E)=E 



(3.2) 



49 



The secular polynomials can be calculated from ^ by taking its kth 
derivative in ^ 

The secular polynomials and energy levels are approximated in this 
chapter using the WKB approximation with initial value conditions 
instead of boundary value conditions (Appendix C). The solutions are 
determined by two free parameters 

WKB M^) 



+ 



<lE - V{x) 

B{E) 
<lE - V{x) 



exp 



exp 



1V2 £ - V{y) dy 
-tV2£ ^E-V{y)dy 



(3.4) 



Initial value conditions have the form ^{xq) — a, -^^{xq) ~ b. For the 
examples discussed in this chapter the generating function has a well 
defined parity and therefore the choice of 



^{0,E)^1 
jlj{-x,E)^tl;{x,E) 



^{0,E)^1 
i/j'iO^E) = 



(3.5) 



as the initial value conditions is obvious. The WKB approximation of 
the generating function for even energy levels is 



WKB 



\ 



E 



E- 



cos 



(3.6) 



The approximate generating function is now used instead of the exact 
function in equation (3.3) to give the approximate secular polynomials 



Pk{E) = 



d'' 



i^^^^i^ix)) 



k\Ak de 



(3.7) 



The results of this approximation were not good enough so the com- 
putation was repeated with higher order WKB approximations. The 
results improved with the order of WKB but the complexity of the 
calculation increased so that large j calculation was still out of reach. 
The results of these approximations are presented in table (3.1) and in 
figures (3.1-3.6) 
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Figure 3.1: The results of differentiation of the harmonic oscillator first 
order WKB approximation j — \- 
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Figure 3.2: The results of differentiation of the sextic potential first 
order WKB approximation j — \- 
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Figure 3.3: The normalized highest energy of the harmonic oscillator 
vs. the order of the WKB approximation, j — \- 
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Figure 3.4: The normalized third highest energy of the harmonic oscil- 
lator vs. the order of the WKB approximation, j = |. 
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Figure 3.5: The normalized highest energy of the sextic potential vs. 
the order of the WKB approximation, j — \- 
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Figure 3.6: The normalized second highest energy of the sextic potential 
vs. the order of the WKB approximation, j — \- 
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Potential 


Exact 


r irst L^rQer vv ivr> 


ThYmit-^-Vi C^rAar 

rourtn ijrcier vvivr) 




14.5 


14.493 


14.5 


Harmonic 


12.5 


12.482 


12.5001 


Oscillator 


10.5 


10.487 


10.5009 




8.5 


8.533 


8.5004 




6.5 


6.639 


6.471 




27.804 


26.711 ±i4.86 


29.168 


Sextic 


17.777 


15.134 ± ill. 046 


23.94 ±i7.187 


Potential 


9.041 


6.455 ± il2.927 


16.012 ±il3.78 




2.273 


±il2.866 


7.85 ± il6.39 




-2.273 


±i5.684 


5.398 



Table 3.1: The exact and approximate energy levels of the harmonic 
oscillator and sextic potentials, in first and fourth order WKB approx- 
imation. 



3.2 Saddle Point Approximation of The 
Secular Polynomials 

The calculation of energy levels through equation (3.7) requires very 
long computations. The fact that the interesting limit is j ^ 1 can be 
used to apply the saddle point approximation on the Cauchy integral 

which is equivalent to equation (3.7). For the saddle point approxi- 
mation of the Cauchy integrals the variables have to be changed into 
new variables that are scaled with j, so that the exponential phase is 
proportional to j. This calculation has many details that are difficult 
to generalize, so two examples are given below. 

3.2.1 Secular Polynomials of the Harmonic Oscil- 
lator 

The harmonic oscillator is chosen to demonstrate this approximation 
because all the energy levels are known exactly. In section (2.3.1) it is 
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shown that for the potential 



V{x) = 

the initial value solution of the Schrodinger equation is 



where 



P,{E)J-\i(E-2i-\ 



Ik 



(3.9) 



(3.10) 



(3.11) 



are the secular polynomials. The zeros of the polynomial P2j+\{E) give 
the first 2j + 1 even energies of the system 



£;r" = 2n+-, n = 0,l,2,.. 



(3.12) 



The approximate WKB initial value solution (3.4) for the harmonic 
oscillator is 



WKB 



E 



cos 



Inserting this definition into the Cauchy integral (3.8) gives 



(3.13) 



Pk{E) « / 

Jc, 



dx e- J 2E 
CO 2m \l 2E-x'^ 



cos 



£ ^2E - 7/2 dy 



(3.14) 



We are interested in the 2j + 1 polynomial. For the proceeding calcu- 
lations it is more convenient to work with the exponential form of the 
cosine 



dx e ^ 



2E 



^1 /co 2m x^^+^ \l 2E-x^ 



(3.15) 



In order to do a saddle point approximation of this integral we need 
a single large parameter (j) that multiplies the entire phase. To this 
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end we will change the variables to scaled variables. The choice of the 
scaling of the variables is such that all the terms in the exponential 
phase will be hnearly dependent on j. The scaled variables are 



y ^ \ljt 

E ^ je (3.16) 
Substituting the scaled variables in (3.15) results in 

^'...WocE^£^/=?;e^-" (3.17) 
where the phase Sa{s) is 



Sais) ^ aiV2 r \l e rft + — -41ns. (3.18) 

JO V 2 2 

This integral will be the starting point of the saddle point approxima- 
tion. The dominant contribution to the integral (3.17) is from constant 
phase paths that pass through the saddle points [31]. The phase Sa{s) 
is expanded to a second order Taylor series around the saddle points s^j 

Sa{s) = S^isp + z)= Sa{sp) + S',isfs)z + \s';^isf,)z\ (3.19) 
The saddle points satisfy 



= m\/2 J e - y + s - ^ = 0, (3.20) 



and are equal to 



s/3(e)=/5jj3^ , /5 = ±1. (3.21) 



The values of the phase and its derivatives at the saddle points are 



Slis,) - -|^. (3.22) 
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We arc interested in paths of steepest descents so we want ^S'^{sg)z^ 
to be real and negative. Since S'^{si3) is negative, z is chosen to be 
reaL The dominant contribution to the integral (3.17) is therefore the 




Figure 3.7: Two paths of integration in the s-plane for the saddle point 
approximation of the harmonic oscillator secular polynomials. 



steepest descent integral 



a,/3=±l 



oJSc 



=2 



2 ^ 



f 

Jo 



(4-£)l 9 



(3.23) 



After integrating and substituting from (3.21) the approximate poly- 
nomial is 



P2,+i(e)« E (4-e)2^ e 

a,/3=±l 



. (3.24) 



Using the parity of the potential to sum over a and /3 the final result 
of the saddle point approximation of the secular polynomial is 



P2,+i(6)oc(4-e)2^ei^^e(4-6) 



cos 



V2e - ^2 dt 



(3.25) 



The real and imaginary parts of F2i+i(e) arc plotted in figures (3.8) 
and (3.9) respectively. The approximate energy levels are given by 
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the zeros of the approximate secular polynomial. P2j+i(e) vanishes for 
e = 0, which corresponds to the lowest exact energy Eq'"'^"' — |. The 
hmit of (3.25) for e — > 4 diverges. The cosine has a zero for e = 2 and 
other zeros, all of which solve the equation 

j J^^W2e - dt ^ {k + ^^ TT. (3.26) 
Expanding the integral gives 

j/^' V2^dt = TT + |(6 - 2) - i(e - 2)=^ - ^ (6 - 2)^ + . . . (3.27) 

The real part of this number is the integral over the allowed region 
j/^ V2e — dt. Solving equation (3.26) in the highest order in 4 
results in 

k = j (3.28) 

This result is valid only to first order (because the previous approx- 
imation was in leading order only), which means that the final result 
of the approximation is 

= 2k (3.29) 

which for large k is an excellent approximation. This result can be 
improved by computing the approximations to higher orders, both in 
the WKB approximation and in the saddle point approximation. 

3.2.2 The Bender-Dunne Polynomials 

The calculation of section (3.2.1) are repeated for the sextic potential 
(1.22), which is the simplest most popular example of QES potentials 

V{x) = ^x''-^{8j + 3)x' (3.30) 
The generating function of the secular polynomials of this potential is 

* = E T^P^i^y (3-31) 

k=Q y^'^)- 
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Figure 3.8: Real part of saddle point approximated harmonic oscillator 
secular polynomial. 




Figure 3.9: Imaginary part of saddle point approximated harmonic 
oscillator secular polynomial. 
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where the secular polynomials Pk{E) are given by the recursion relation 

Pk+i{E) = EPk{E) - 2k{2k - l)(2j + 1 - k)Pk-i{E) (3.32) 

with the initial conditions Pq = 1, Pi = E. The first energy levels and 
secular polynomials are listed in table (1.1). The WKB approximated 
generating function is 



jWKB _ J E 



V2 ^Je - V{y) dy 



(3.33) 



Incorporating this into (3.8), taking k = 2j + 1 and writing the cosine 
as a sum of integrals gives 



^4 



The contour integrals can be approximated using the saddle points 
method for large j. The saddle point approximation requires us to de- 
fine scaled variables so that the exponential phase is linearly dependent 
on j 

1 

X = j-^s 

y = j'^t 

E = re 

V{x) = -l(8j + 3)^2 = ji[i/- 4^2]. =ji^;(s) (3.35) 
Substituting these variables in P2j+i{E) leads to 



P. 



where the phase is 

= aiV2 r Je-vit) dt+^— -Alns. (3.37) 
Jo 4 
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This is the starting point of the saddle point approximation. To find 
the integration paths the phase Sa{s) is now expanded in a Taylor series 
around the saddle points s^j 

Sa{s) = S^{sf,{e) + z)^ S^is^ie)) + ^S';,{sp{e))z' . (3.38) 
The saddle points are the points where the first derivative vanishes 

-^Sa = aiV2^e - v{s) +s^-- = Q. (3.39) 
The saddle points of the integral (3.34) are 

s^(e)=/3z^, /3 = ±1. (3.40) 

The values of 5'^ and its derivatives at the saddle points are 

fsi3ie) I 16 64 

^a(s/3(e)) = ^^^^/o \Je-v{t)dt + — -\n{—){;iAl) 

^/a{s^{e)) = (3.42) 

^.S.(s,ie)) = (3.43) 

The range of epsilon for the QES sector is 

16/2 16/2 

-yV3<^<yV3' ^'-''^ 

The second derivative must be real and negative for paths of steepest 
descents. Since the second derivative changes sign at e = ±4 and the 
interesting energies are the highest, the range is limited to 

4<e<yy|:^ 4.35 (3.45) 

though the final result should be the same for the entire range of e. 
For the chosen range the path parameter z has to be real and the 
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integration paths are plotted in (3.7). The main contribution to the 
value of P2j+i (e) is therefore from the integrals 



4 POO 



dz e ■'■^(e-ie)^ . (3.46) 



After the integration, of (3.40) is substituted and the parity of the 
potential is used to get the approximate value of P2j+i(e) 



2i+l 



^e) oc e 



^^+'e^'^ cos 



iV2 



2V2 



v{t) dt 



(3.47) 



The approximate polynomial can be evaluated numerically as shown 
in figures (3.10-3.11). The zeros of this function are the zeroes of the 
cosine, which can be calculated order by order in j. The leading order 
calculation results in 



16 /2 .3 
= — \ - 72. 
3 V 3 



(3.48) 



which is the approximate highest energy of the QES sector obtained in 
[21]. The numerical results are summed in table (3.2). This is a first 
order approximation of the highest QES energy that is obtained through 
a general calculation independent of the energy reflection symmetry. 



J 


Exact 


Saddle Point Approximation 


15 
2 


88.404 


89.443 


8 


97.461 


98.534 


17 

2 


106.809 


107.915 


9 


116.438 


117.576 


19 
2 


126.34 


127.508 


10 


136.507 


137.706 



Table 3.2: The exact and saddle point approximate energy levels of the 
harmonic oscillator. 
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1x10 
7.5x10 




-1x10 



Figure 3.10: Real part of saddle point approximated Bender-Dunne 
polynomial. J = |- 




Figure 3.11: Imaginary part of saddle point approximated Bender- 
Dunne polynomial. J = |- 
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3.3 Summary of Results 



The WKB approximation of the generating function is a new approach 
to the approximation of energy levels. Combined with the saddle point 
approximation the complexity of the energy calculation is small. This 
calculation is a first order approximation of the energy levels of the QES 
, potentials to leading order in \. For the sextic potential the calculation 
approximates the highest QES sector energy level as a function of j 
without any dependence on the energy reflection symmetry. 
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Chapter 4 

Summary and Conclusions 



Quasiexactly solvable potentials have several unexplored properties. 
This work concentrates on the connection between initial value and 
boundary value solutions of the one-dimensional Schrodinger equation, 
which is then used as a basis for approximation. 

Quasigauge Hamiltonians in one dimension can be constructed from 
any bihnear combination of the generators of the sl{2) Lie group. The 
basis of the solutions is constructed from the jth representation of the 
group and an arbitrary set that spans the orthogonal space. As a result 
the quasigauge Hamiltonian has a block diagonal form with at least 
one finite block which has the same dimension as the representation 
{2j + l). The wave function that solves the Schrodinger equation has an 
exponential factor and a pre-exponential function which is a polynomial 
of degree 2j + 1 in the group variable ^. The Schrodinger equation may 
be replaced by a recursion relation for the wave function coefficients. 

In chapter (2) exact solutions of some QES potentials and some 
exactly solvable potentials are discussed. In section (2.1) a new so- 
lution of the most general one-dimensional quasigauge Hamiltonian is 
constructed. This solution includes the transformation between the 
Hamiltonian and the quasigauge Hamiltonian Hq by calculating the 
change of variables, the exponential behavior of the wave functions and 
the potential. The recursion relations for the wave function coefficients 
is written in matrix form. The quantization condition for this system 
is given by the vanishing of the secular polynomial, which is the de- 
terminant of the coefficient matrix. The secular polynomials are found 
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to be related to each other, through a general recursion relation, with 
initial conditions. This simplifies the computation of the energy levels 
for any one-dimensional QES potential. The coefficients of the pre- 
exponential part of the wave function are found to be proportional to 
the minors of the coefficient matrix, which explains the connection be- 
tween the boundary value conditions and the quantization condition. 
The connection between the even and odd parity potentials of the same 
family (i.e. with the same i{x)) is shown. In section (2.2) it is proved 
that for all QES potentials with €k — 0, the generating function of 
the secular polynomials is the initial value solution of the Schrodinger 
equation, thus expanding the results of existing works. In section (2.3) 
the one-dimensional exactly solvable potentials are analyzed using QES 
formalism and the difference between exact solvability and quasiexactly 
solvability is discussed. In section (2.5) new recursion relations for the 
secular polynomials of the SUSY potential are derived. The separate 
recursion relations for coefficients of the upper and lower parts of the 
spinor wave function are found. 

In chapter (3) the exact energy levels of QES potentials are approx- 
imated using a new technique in order to simplify the calculation for 
large j. The generating function of the secular polynomials is approxi- 
mated using the WKB approximation with initial value conditions (in 
contrast to the more common boundary value WKB solution). The 
approximate energy levels are derived by replacing the exact generat- 
ing function by the approximate function. Since the complexity of this 
calculation is considerable, a second approximation is introduced. The 
derivatives of the generating functions are replaced by Cauchy integrals, 
which arc then calculated using the saddle point approximation. This 
method is demonstrated for two specific potentials - the harmonic os- 
cillator and the sextic potential. The zeros of the approximate secular 
polynomials (which are the approximate energy levels) are calculated 
order by order in the group parameter 4. For the harmonic oscillator, 
for an energy of order j the error is of order 1. The results arc valid for 
energies that are not in the top part of the calculated spectrum, because 
the approximate secular polynomial has singular points for high energy 
levels. For the sextic potential only the highest energy of the QES part 
of the spectrum can be evaluated. For the highest energy of order j^^^ 
the error is of order j^^^. This result was obtained before by numeri- 
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cal calculations and proved analytically based on the energy reflection 
symmetry. Using our approximation this result can be duplicated for 
any potential of interest without this symmetry. 

This work can be continued by generalizing the calculations to in- 
clude the higher dimensions and other groups, including SUSY groups. 
The approximation can be applied to more potentials and improved 
by introducing the saddle point approximation to higher order, with 
higher order WKB approximation. 
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Appendix A 
Oscillation Theory 



In one-dimensional potentials the energy levels are not degenerate and 
therefore can be ordered by magnitude. The wave functions can be 
numbered according to the number of zeros (or nodes). The two num- 
bering methods coincide. 

In the solution of the one-dimensional eigenvalue problem 

H{x)i){x) = -tlj"{x) + V{x)^{x) = Eip{x) (A.l) 

there is no degeneracy of the energy levels, since if two different solu- 
tions ipi and ip2 

-ij'l + Vij, = Eii,^ (A.2) 
-ijl + Vij2 = £^2V'2 (A.3) 

are assumed to have the same energy 

= E2, (A.4) 

then the difference ~ be equal to zero 

-<V2 + V"^1^2 = ^^1^2 (A.5) 

-V^^Vi + l^V'aV'i = Ei}2iJi (A.6) 

- (V'lVs - V'iV'2) = -^(V;V'2-V'2V'i) = 0, (A.7) 
which means that the Wronskian W{x) 

W{x) = ip[ip2 - '02'0i = const (A. 8) 
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is constant. The Wronskian must vanish for x ±00 so it is constantly 
zero 

W{x) ^^~= ^ W{x) = Q ^ Ei^ E2. (A.9) 

This means that the two wave functions must be identical or the ener- 
gies not equal. 

Repeating the calculation (A.2-A.7), with E2 > Ei results in 

- = (E, - E2)i^llP2. (A.IO) 

If ■^i has zeros in a; = a, 6 and is positive in the section (a, b) then 
integrating (A.IO) from a to 6 leads to 

- (V'iV2 - i^2i^l)t = (^1 - E2) /%iV2 dx. (A.ll) 

J a 

Since ipi{a) = 0, 'ipi{b) = 0, ipi^a < x < b) > 0, the derivatives of ipi 
at the end points are known to be ip'i{a) > 0, ip'i{b) < 0. If we assume 
ip2{(i) > then the left hand side of (A.ll) is positive 

jC = iP[{a)Mo') - ip'iib)Mb) > 0, (A.12) 

and if we also assume that ip2{(i < x < b) > then the right hand side 
is negative 

n={Ei- E2) [ dx < 0, (A.13) 

which is a contradiction. The assumption that ip2{0' ^ x < b) > is 
proved to be wrong and therefore ip2 has a zero in the section (a, 6). 
Since at the boundary both wave functions vanish, ip2 must have one 
zero more than ip2- Since the ground state has no zeros, the number of 
nodes of the n-th eigenfunction (corresponding to the n-th eigenvalue) 
is n. A more thorough discussion on oscillation theory can be found for 
example in [32, 33]. 
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Appendix B 

Initial Value and Boundary 
Value 

A linear differential equation has a unique initial value solution and any 
number of boundary value solutions. In an eigenvalue problem only for 
certain values of the free parameter the system has boundary value 
solutions. For distinction between the initial value and boundary value 
solutions of the same equation ^ is used for the initial value solution 
and ip for the boundary value solutions. 
The equation 

{-d^ + W)^ = (B.l) 
has a single solution "i/x that satisfies the initial conditions 

^x{t = 0) = 0, dt^xit = 0) = 1. (B.2) 

On the other hand, the operator —df + W that acts on the space of 
functions that obey the boundary conditions ijj{t — 0) — 0, ijj{t — T) — 
has an eigenvalue A„ if and only if 

V^A^W^O, V'A^m^O. (B.3) 

The determinant of the operator —d^+W is defined to be 

det{-d^ + W)^l[Xr,. (B.4) 

n 
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It is proved in [34] that for two different functions H/^^'^) witli initial 
value solutions ' the equation 



det 



-df + - A 



(T) 



(B.5) 



is satisfied. As a result the normalization factor defined by 

det(-a2 + W) 



*A=o(r) 



(B.6) 



is constant. Applying this theorem to nth order difference equations 
[31], the equation 

MfcxfePfe = APfe (B.7) 

— * 

has a single solution Pfe(A) that satisfies the initial conditions Pm — 

0, Pm+i = 1- On the other hand, the matrix Mkxk that acts on C^, 
which obeys the boundary conditions Cm = 0, Cm' = for m — m' ^ 
0, ±1, has an eigenvalue A„ if and only if Cm(An) = 0, Cm' (An) = 0. 
The determinant of the matrix M^xk is defined to be 



detMjfcxit = n-^"- 



(B.8) 



For two different matrices M^^fJ with initial value solutions Pk the 
equation 

^iV(A) 



det 



(B.9) 



with Ik the k X k unit matrix, is satisfied. This theorem can be proved 
following the proof given in [34] for the differential operators. As a 
result the normalization factor A^, defined by 



det Mkxk 



(B.IO) 



Pm'{X = 0) 

is constant, and the zeros of Pm'{X — 0) are the eigenvalues of det M^xk- 
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Appendix C 

WKB Approximation 



The WKB approximation is a semi-classical approximation of the wave 
functions and energy levels of the Schrodinger equation. In this work 
the well knows WKB approximation [4, 31, 35] is expanded to include 
initial value solutions of the Schrodinger equation. In the WKB ap- 
proximation the solution of the Schrodinger equation 



for small h is assumed to have the general form 

\WKB 



exp 



-| oo 
"' n=0 



(C.l) 



(C.2) 



Since the Schrodinger equation is a second order equation the general 
solution has two free parameters. 

^WKB ^ j^(^^^^i{n-^So-ns2+n^S4—-)^(Si-n'^S3+n'^S5—-) 

The terms Sn are found by solving equation (C.l) order by order in h. 
The exphcit solution to first order in h is 

A{E) 



WKB 



+ 



E - V{x) 
B{E) 
P - V{x) 



exp 



exp 



iV2£p-V{y) dy 
-iV2j^^jE-V{y)dy 



(C.4) 
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For any order of the approximation two more conditions are required 
to determine A{E) and B[E). The quantization of the wave functions 
and energy levels is a result of the boundary value conditions, which 
are given by the normalization condition 

ip{x ±00) ^ 0. (C.5) 

The exact and approximate initial value solutions of the Schrodinger 
equations have initial conditions of the form 

^>(x = 0) = a, ^-^ix = 0) = b. (C.6) 
ax 

If the potential is even then the generating function will have a well 
defined parity. Therefore the initial conditions for even potentials will 
be either 

^(x = 0) = 1, = 0) = (C.7) 

ax 

or 

^(x = 0) = 0, ^"^{x = 0) = 1 (C.8) 

depending on the parity of the states. The result of this approach is 
the initial value WKB approximation. 
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